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Abstract 

The problem of existence of non-analytic (Griffith-like) contribu- 
tions to the free energy of weakly disordered Ising ferromagnet is stud- 
ied from the point of view of the replica theory. The consideration 
is done in terms of the usual random temperature Ginzburg-Landau 
Hamiltonian in space dimensions D < 4 in the zero external magnetic 
field. It is shown that in the paramagnetic phase, at temperatures 
not too close to Tc (where the behaviour of the pure system is cor- 
rectly described by the Gaussian approximation), the free energy of 
the system has additional non-perturbative contribution of the form 
exp{— (const) ^!^—-^^} (where t = {T — Tc)/Tc), which has essential 
singularity in the parameter u ^ which describes the strength of 
the disorder. It is demonstrated that this contribution appears due 
to non-linear localized (instanton-like) solutions of the mean-held sta- 
tionary equations which are characterized by the special type of the 
replica symmetry breaking . It is argued that physically these replica 
instantons describe the contribution from rare spatial "ferromagnetic 
islands" in which local (random) temperature is below Tc. 
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1 Introduction 



According to the original statement of Griffith [|I[] , the free energy of the ran- 
dom Ising ferromagnet in the temperature interval above its ferromagnetic 
phase transition point Tc and below the critical point TjP'' of the correspond- 
ing pure system must be a non-analytic function of the external magnetic 
field h, such that in the limit — > the free energy as a function of h 
have essential singularity. Since this type of phenomena, namely, existence of 
non-analytic non-perturbative contributions to thermodynamical functions 
in random systems, seems to be rather general one, at present it became 
common to call any such contribution as the " Griffith singularity" . 

Due to intensive theoretical [@] and numerical studies of the Grif- 
ffih singularities it was also discovered that the dynamical properties of the 
system in the temperature interval <T < Tj°) are not just ordinary para- 
magnetic. In particular, the time correlation functions here are described 
by the so-called stretched-exponential asymptotic behaviour which is much 
slower than the usual exponential one, as it should be in the paramagnetic 
phase. To underline that the properties of the system in the temperature 
interval < T < T]°) are not quite paramagnetic, it became common to 
call the state of the system here as the " Griffith phase" . 

At the level of "hand- waving arguments" the dynamical Griffith phe- 
nomena can be explained "theoretically" rather easily: considering e.g. the 
bond diluted Ising model, one can note that at temperatures below Tj^-* in 
the "ocean" of the zero magnetization paramagnetic background the random 
system must contain disconnected locally ordered "ferromagnetic islands" 
(composed only of the pure system bonds) of all sizes, which, in turn, creates 
the whole spectrum (up to infinity) of relaxation times. Having an infinite 
spectrum of relaxation times, with some imagination, it is not difficult to 
derive any relaxation law one likes, and the stretched-exponential one in 
particular. 

Although it is commonly believed that the main point of the above "ex- 
planation" , namely the existence of infinite number of local minima states, 
must be a general key for understanding the Griffith phenomena (both dy- 
namical and statistical mechanical), despite many efforts during last thirty 
years, it turned out to be extremely difficult to construct more or less elabo- 
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rated and convincing theory. For that reason any progress in understanding 
of the effects produced by numerous local minima states (which, so to say, 
are away from the perturbative region) looks valuable. 

In this paper I am going to study non-perturbative contributions to the 
thermodynamical functions of weakly disordered (random temperature) D- 
dimensional {D < 4) Ising ferromagnet in the paramagnetic phase away from 
the critical point. In the continuous limit this system can be described by 
the usual Ginsburg-Landau Hamiltonian: 



H[(t){x)]5T{x)] = J d^x 



\ (V0(x))^ 



5t{x)) <P^{x) + -g<p\x) 



(1.1) 

Here r = (T — Tc)lTc <^ 1 is the reduced temperature, and the quenched 
disorder is described by random spatial fluctuations of the local transition 
temperature 5t{x) whose probability distribution is taken to be symmetric 
and Gaussian: 



P[5r] =poexp(^-i-|ci^x(5r(a;))^j , (1.2) 

where u <C is the small parameter which describes the strength of the dis- 
order, and po is irrelevant normalization constant. For notational simplicity, I 
define the sign of 5t{x) in eq.(p,.l|) so that positive fluctuations lead to locally 
ordered regions, whose effects will be the object of our further study. 

As far as the corresponding pure system {u = 0) is concerned, it is well 
known that in the close vicinity of Tc, at |r| <^ ~ g'^/^'^~^\ its properties 
are defined by non-Gaussian critical fluctuations (which can be studied e.g. in 
terms of the e-expansion renormalization group approach), while away from 
Tc, at |r| ^ Tg, the situation is getting Gaussian, and everything becomes 
very simple. Here the total magnetization of the system is defined by the 
order parameter (0) = (poij) which is equal to above Tc, and it is equal 

to ±y^jrj/^ below T^, the asymptotic behaviour of the correlation function 
G{x — x') = ((0(x)0(x')) — 0o) is defined only by the Gaussian fluctuations: 
G{x) ~ and the singular part of the free energy /(r) scales with 

the temperature as /(r) r^/^. 



3 



Usually, the random system, defined by the Hamiltonian ( p. .11 ) , were stud- 
ied from the point of view of the effects produced by the quenched disorder 
on the critical phenomena in the close vicinity of the phase transition point. 
Renormalization group consideration shows that if the temperature is not 
too close to Tc, at Tu <^ T <^ Tg (where the disorder dependent crossover 
temperature scale ~ m^^" is defined by the specific heat critical exponent 
a > of the pure system) the critical behaviour is essentially controlled 
by the pure system fixed point, and the disorder produces only irrelevant 
corrections. On the other hand, in the close vicinity of the critical point, at 
T -C Tu, the critical behaviour moves into a new universality class defined 
by the so called random fixed point, which turns out to be universal [^. In 
recent years, however, this very nice physical picture has been questioned on 
the grounds that the renormalization group approach completely misses the 
presense of numerous local minima configurations of the random Hamilto- 
nian ( p. . 1| ) , which, in principle, may cause the spontaneous replica symmetry 
breaking in the interaction parameters of the critical fluctuations, which, in 
turn, may ruin the above physical scenario 0. 

Leaving the discussion of this very difficult problem for future analycis, 
in this paper I would like to pose much more simple question: how the ther- 
modynamic functions of this system depend on the strength of the disorder 
u (in the limit u 0) far way from at r ^ Tg, where the behaviour of the 
pure system is correctly described by the Gaussian approximation? It tuns 
out that even this, as if almost trivial question is not so easy to answer. 

Of course, first of all, one can proceed in a straightforward way, developing 
the perturbation theory in powers of the parameter u at the background of the 
pure system paramagnetic state (0) = using the Gaussian approximation 
for the thermal fluctuations. There is nothing wrong in this approach, but the 
problem is that it can not give all thermodynamic contributions which exist 
at M 7^ 0. The drawback of this type of the perturbation theory is the same 
as that of the renormalization group: it completely misses the existence of 
numerous (macroscopic number) local minima configurations of the random 
Hamiltonian ( |1 . 1| ) . 

At the level of "hand waving arguments" it is very easy to see what all 
these off-perturbative states are. At any u ^ there exists finite (expo- 
nentially small) density of "ferromagnetic islands" in which local (random) 
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temperature is below Tc (such that St{x) > r), and the minimum energy 
configurations here are achieved at non-zero local value of the order parame- 
ter: (poix) ~ ±^ {6t — T)/g. Since the spatial density of such islands is finite, 
and each island provide two (±) possibilities for the local magnetization, 
the total number of the local minima configurations in the system must be 
exponential in its volume. 

Formally, to take into account the contributions of all these states, one 
has to proceed as follows. For an arbitrary quenched function St{x) one has 
to find all possible local minima solutions of the saddle point equation: 



A0(x) + (r - (5r(x))0(x) + (j){xf = (1.3) 



Then one has to substitute these solutions into the Hamiltonian ( |1.1| ) and 
calculate the corresponding thermodynamic weights. Next, to compute the 
partition function one has to perform summation over all the solutions, and 
finally to get the corresponding free energy one has to take the logarithm 
of the partition function and average it over random functions 6t{x) with 
the probability distribution (|1.2|) . Clearly, it is hardly possible that such a 
programme can be implemented. 

On the other hand, as usual, for the systems which contain quenched 
disorder we can use the standard rephca method and reduce the problem 
of the quenched averaging to the annealed one for n copies of the original 
system: 

F = -(hrZ) = -limi(z^-l) (1.4) 



where (...) denotes the averaging over random functions St(x) with the prob- 
ability distribution (|1.2|) , and 

Z[6t{x)] = J V(j){x) exp {-H [0(x); r(x)]) (1.5) 

is the disorder dependent partition function which is given by the functional 
integration over configurations of the field (f){x). 

Simple Gaussian integration over 6t{x) in eq.( |1.4|) yields: 
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n 

a=l 



exp {-Hn[(f)aix)]) 



;i.6) 



where 



1.7) 

is the spatially homogeneous rephca Hamiltonian. 

Now, if we are intended to take into account non-trivial local minima 
states, instead of solving the original inhomogeneous stationary equation 
(|1.3|), we can consider the corresponding replica saddle-point equations: 



- A(j>a{x) + tMx) + ^l{x) - uMx) E = (1-8) 

6=1 

Since until now all the transformations were exact, these equations must 
contain (may be in a slightly hidden way) all the relevant non-trivial states 
which in the language of the original random Hamiltonian correspond to rare 
ferromagnetic islands. 

At this stage we can note one very simple point. Looking for various types 
of solutions of the above equations one can try first of all the simplest possible 
"replica symmetric" ansatz, in which the fields in all replicas are assumed 
to be equal: (f)a{x) = <P{x). In this case the last term in the eqs.( |1.8|) (which 
contains the factor Ylb=i (Pbi^) = n(j)'^{x)) drops away in the hmit n — >• 0, and 
these equations reduce to the pure system saddle-point equation: 



- A(j){x) + T(j){x) + (j){xf = (1.9) 

which at r > has only trivial solution 0(x) = 0. It means that in any 
non-trivial solution of the eqs. (|1.8|) the fields (f)a{x) in different replicas can 
not be all equal. In other words, the symmetry among replicas in the replica 
vector (f)a{x) must be broken. 
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The methodological aspects of how to handle with the vector replica sym- 
metry breaking situation in various disordered systems are described in the 
paper In the next Section this method will be applied for the problem 
described above. It will be shown that indeed, in the high-temperature re- 
gion (r > 0) eqs. (|1.8| ) have non-trivial localized (having finite size and finite 
energy) solutions in which the replica symmetry in the fields (pai^) is broken. 
The formal summation over all such solutions provides the contribution to 
the free energy of the typical Griffith- like form: exp{~{const)^^^—^^}. It will 
also be shown that the mean-field approach (in which the critical fluctuations 
are ignored) used in this paper is grounded only if the temperature is not 
to close to Tc, namely at r S> ~ g'^^^'^~^\ the same as in the classical 
Ginsburg-Landau theory. Finally, it will be demonstrated how this type of 
non-analytic contribution to the free energy can be estimated from purely 
physical arguments taking into account probabilities for the typical "ferro- 
magnetic islands". 

To avoid possible misunderstandings, in conclusion of this introductory 
Section I would like to note the following essential point. The problem con- 
sidered in this paper is actually rather far from the original one studied by 
Griffith as well as by many other people later on. Since the shift of Tc in the 
weakly disordered ferromagnet compared to Tj°) of the pure system is of the 
order of ^/u, in the limit u <^ g the interval of temperatures Tc < T < Tj°^ 
where the so called Griffith phase is expected to take place, appears to be 
well inside of the temperature interval Tg ~ (^2/(4-d) ^j^g^-g ^j^g critical fluc- 
tuations are essential, and where the mean field approach considered in this 
paper can not be used. For that reason, in the considered range of temper- 
atures T ^ Tg it is hardly reasonable to look for non-analytic behaviour of 
the free energy as the function of the external magnetic field (at least the 
present approach in terms of the replica instantons modified by the external 
field h does not seem to indicate on any non-analyticity in h). The aim of 
this paper is just to demonstrate that in additional to the "usual" Griffith 
singularities in terms of the external field, the free energy of the random fer- 
romagnet (in the zero magnetic field) must also be non-analytic in the value 
of the parameter which describes the strength of the disorder. 
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2 Replica Instantons 



Following the general strategy developed in the the paper 0, let us assume 
that in addition to the trivial replica symmetric (RS) solutions of the saddle- 
point equations ( p..8| ) there exist other types of solutions, which are well 
separated in the configurational space from the RS state. In this case, denoting 
the contribution of these non-trivial states by the lable "replica symmetry 
breaking" (RSB), the replica partition function, eq.(|L6|), can be decomposed 
into two parts: 

Z^ = Zns + ZnsB (2.1) 

where Zj^s contains all the perturbative contributions in the vicinity of the 
trivial state (f)a{x) = 0. As usual, this partition function can eventually be 
represented in the form: 

Zrs = exp i-nVfRs) (2.2) 

where V is the volume of the system, and fjis is the free energy density, 
which contains the pure system leading term ~ r^/^ (at temperatures not 
too close to Tc, r 3> Tg), plus the perturbation series in powers of the disorder 
parameter u. 

Thus, in terms of the general replica approach, according to eq. 
the total free energy we get: 

F = Vhs + Frsb (2.3) 
where the additional RSB part of the free energy 

Frsb = - lim -Zrsb (2.4) 

must contain all non-perturbative contributions (if any) which are away from 
the trivial state (pa = 0. It is this part of the free energy which will be point 
of our further study. 

The simplest possible non-trivial replica structure for the solutions of the 
saddle-point equations ( |1.8D can be taken in the following form (see |0): 
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(2.5) 



for a = 1, ...,k 
foT a = k + 1, ...,n 

where k is the integer value parameter: k = l,2,...,n which defines a given 
structure of the trial replica vector (pa (note that the value k = should be 
excluded since it describes the trivial RS solution which is already taken into 
account in /rs)- 

Substituting this anzatz into eqs. (|1.8|) as well as into the replica Hamilto- 
nian ( p..?] ), one finds that for a given value of the parameter k the fields (j){x) 
in eq. (|2.5| ) are defined by the solutions of the following saddle-point equation: 



- A(j){x) + T(t){x) - X{k)(j){xf = (2.6) 
and the thermodynamic weight of any such solution is defined by the energy: 



E{k) = k j d^x 



\{yct>{x)f + \T<p\x)-\m(p\x) 



(2.7) 



where 



\{k) = {uk-g) (2.8) 

Summing over the parameter k and taking into account the combinatoric 
factor which is the number of permutations among replicas in the ansatz 
structure (P75|) for the free energy, eq. (^.41), one gets: 



1 " n' 

Frsb = - lim - 5: -— — — exp{-E(fc)} (2.9) 
n kl[n — k)\ 

To take the limit n — >■ the series in the above equation can be be represented 
as follows: 



1 °o r(n 4- 1) 

Frsb = - lim - ^ ^ exp{-i?(A:)} (2.10) 

n ^ L [k + 1)1 [n — k + 1) 

Here the summation over k is extended beyond k = n to oo since the gamma 
function is equal to infinity at negative integers. Now using the relation 
T{—z) = 7c[zr{z) sin(7rz)]~-^, we can perform the analytic continuation n — > 0: 
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r(n+l) ^ 

r(fe+i)r{n-fc+i) 

_ r(n+l)(fc-l-n)r(fc-l-ra) sin(7r(fc~l-n)) i ^ (^l)''-^ 

~ 7rr(/c+l) l(n-»0) — k 

Thus, for the free energy ( p.9|) one obtains: 



(2.11) 



Frsb = - E exp{-E{k)} (2.12) 

fc=i 

At this stage we can note the following important point. For any non 
localized (e.g. space-independent) solution, such that its energy ( |2.7| ) is di- 
vergent with the volume V of the system, the corresponding contribution 
to the free energy ( p.l2|) will not be proportional to V, but instead, it will 



contain the volume in the exponential factor. It means that at least for the 
bulk properties of the system this type of solutions must be irrelevant. 

Thus, we have to look for localized solutions: the ones which are local 
in space (breaking translation invariance) and which have finite (volume- 
independent) energy. Let us suppose that such instanton-type solutions do 
exists (see below) , and that for a given k the solution is characterized by the 
spatial size R{k). Then, if we take into account only one-instanton contribu- 
tion (or in other words if we consider a gas of non-interacting instantons), 
due to obvious entropy factor V/ (which is the number of positions of the 
object of the size R in the volume V) we get the free energy proportional to 
the volume: 

oo r -| \A;— 1 

Frsb ^-VJ2 ^-^—R-^'ik) exp{-E(A;)} (2.13) 
fc=i 

Now let us come back to the saddle-point equation ( |2.6| ), and let us con- 
sider the range of the parameter k such that \{k) = {uk — g) > (i.e. 
k > [g/u]). Rescaling the fields: 



<j>{x) = ± /^^(xv/7) (2.14) 
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instead of eq.( |2.6| ) one get the following differential equation which contains 
no parameters: 



Correspondingly, for the energy, eq. (|2.7]), one obtains: 



(2.15) 



where 



E{k) 



k 



uk — g 



(2.16) 



^ivi^{z)r+^^'{z)-\ij\z) 



(2.17) 



The equation (|2.15|) is well know in the field theory (see e.g. 0): it is for 
the present choice of signs of the linear and the cubic terms (which imposes 
the conditions: r > and k > [g/u]) in dimensions D < 4 this equation has 
spherically symmetric instanton-like solutions such that: 



ip{\z\ < 1) ~ ^(0) ~ 1, 

(2.18) 

ij{\z\ > 1) ~ exp(-|2;|) 0. 

The energy, eq.( [^.17| ), of such a solution is a finite and positive number. Of 
course, for a generic value of the field ip{0) at the origin, the solution tends 
to the values ipil^l C)o) = ±1 which are the extrema of the potential 
^ip^ — \ip^, and any such solution has divergent with volume energy (|2.17|) . 
However, there exists a discrete set of initial values ipQ for which the solutions 
(exponentially) tends to zero at infinity, and which have finite energies. It 
can be show that the solution with the minimal energy Eq corresponds to 
the minimal value of |^o| in the set. In particular, at D = 3, t/'o — 4.34 and 
Eq ~ 18.90. For our further calculations with the exponential accuracy it will 
be sufficient to take into account only the solution with the minimal energy. 

According to the rescaling (|2.14|) , in terms of the original fields 0(x) the 
size of the instanton is i? = r~^/^ (note that it does not depends on /c), which 
coincides with the usual correlation length of the Ginsburg-Landau theory. 
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Substituting this value of R as well as the energy ( |2.16| ) of the instanton into 
the series (|2.13| ) for the free energy one gets: 



RSB 



oo 

E 

k>[g/u] 



\k-l 



k 



-2'^exp 



uk 



-T 



(A~D)/2 



(2.19) 



(the factor 2^ appears due to independent summation over ± signs, eg . (12.141) , 
in k no n- zero replicas, eq.( |2.5| )). It can be easily shown that under the 
considered conditions on the parameters u,g and t [u <^ g <^ 1, and 
g2/{A-D) <^ 7" <^ 1) the leading contribution in the above series with the 
exponential accuracy comes from the region k ^ g/u ^ 1: 



V 



RSB 



n II 
~ r ' exp 



1-(4-D)/2' 



U 



X E 



:-i 



\k-l 



k 



(2.20) 



Here the absolute value of the series Y^k'=ko>i k~^{—lY~^2'' can be estimated 
by the upper bound ~ k^^2''°, and since it is assumed that r ^ g'^/i^^^) the 
term ^ In 2, which appear in the exponential, can be dropped in comparison 

(4 — 1?) /2 

with Eq- — - — . Thus, for the density of the free energy we finally obtain the 
following contribution: 



V 



RSB 



exp 



~E> 



T 



(4-D)/2 



U 



(2.21) 



(where we drop all pre-exponential factors, which within the present accuracy 
of calculations can not be defined). 



3 Fluctuations 

Note first of all, that one should not be confused by the "wrong" sign of the 
0^ interaction term in the energy function ( p.7|) , which for the usual field 
theory would indicate on its absolute instability. Here, as usual in the replica 
theory, in the limit n ^ everything turns "up down", so that the minima 
of the physical free energy actually correspond to the maxima of the replica 
free energy. It can be easily shown (see below) that formal integration over 
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n-component replica fluctuations around considered instanton solution in the 
limit n — i> yields physically sensible result. 

Proceeding the same way as in the usual Ginsburg-Landau theory, let 
us determine under which conditions the above mean-field approach used to 
derive the result (|2.21| ) can be valid. Introducing small fluctuations ipa{x) near 
the instanton solution, eqs.( |2.5| ) and ( |2.14| ): (paix) = <^^*"**^(a;) + ipa{x), in the 
Gaussian approximation we get the following Hamiltonian for the fluctuating 
flelds: 



H[<f] = J d^x 



a=l 



a,b=l 



where the matrix Tahix) contains the k x k block: 



(3.1) 



1 - ^^^^^^(xv^)l 6ab - -^^\x^) (3.2) 



uk — g J uk — g 

(where a,b = 1, ...,k) and the diagonal elements for the remaining (n — k) 
replicas: 



T. 



[n-k) 



uk 



ab 



-^l^\xV^) 6, 



-r y- M ■ , , -^ab (3.3) 

y uk — g J 

(where a,6 = + 1, ...,n). Here the function il){z) is the instanton solution, 
eq-(|2J[g). 

Since the mass term in the Hamiltonian ( p.l|) is proportional to r, the 
behaviour of the correlation function of the fluctuating flelds at scales 
|x| -C -Rc ~ T"^''^ appears to be the same as in the Ginsburg-Landau theory: 
Gab{x — x') = {Lpa{x)Lpb{x')) ~ |x — x'\^^^^'^'> Sab (bcyoud Rc this correlation 
function decays exponentially). Therefore, the typical value of the fluctua- 
tions {(f"^) can be estimated in the usual way: 



-jZR-^ I d^xGUx)-T^^-'^/' 



(3.4) 
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The saddle-point approximation considered in the previous Section is justified 
only if the typical value of the fiuctuations is small compared to the value of 
the "background" instanton field 0(*"'^*)(a;) ~ ^r/X{k) (see eq. (|2.14|) ): 

A(A;) 



On the other hand, the contribution ( p. 211) appears due to summation in 
the region k ^ g/u. Thus, one can get this type of contribution to the free 
energy only in the following interval of summation over k: 



y « A; « Lr^^-D)/2 (3.6) 
u u 

This interval exists provided 



r » g^/^'-''^ (3.7) 

which is the usual Ginsburg-Landau criteria. 

One can also arrive to the same conclusion deriving the fiuctuational 
contribution to the RSB part of the free energy by the direct integration 
over fiuctuating fields using Gaussian Hamiltonian ( |3.1| ) (this way one can 
also check that this contribution contains no imaginary parts which would 
happen, if the considered extrema would correspond to physically unstable 
field configuration). Assuming the 6'-like structure of the instanton solution: 
'^(kl < 1) — ^(0) = ■?/'o ~ 1 and '0(|;2| > 1) = 0, the fiuctuating modes with 
momenta p <^ ^/t and p ^ i/r in the Hamiltonian (|3.1|) can be explicitly 
decoupled: 



H 



1 " r d^p 



P^^ab + rTab '^aip)<^bi-p) 



2 ^^;^7|p|<v^ {2tt 



d^p 



P'^lVaip)? 



(3.8) 



where p-independent matrix Tab is given by eqs.( |3.2| )-(|3.3|), in which instead 
of the function ip^Xy/r) one has to substitute the constant ipo- 
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The integration over the replica symmetric modes with momenta p <^ i/r 
(they corresponds to fluctuations at scales much bigger than the size of the 
instanton) described by the second term of the Hamiltonian ( |3.8| ), gives the 
contribution of the form ex]i{—nVfRs), and it vanishes in the limit n — 
(note that in the RSB part of the free energy we have to keep only the terms 
which remain finite in the limit tt, — >■ and not linear in n). This is natural, 
because this contribution is already contained in the RS part of the free 
energy. 

The intergation over the modes with momenta p ^ ^/r is slightly cum- 
bersome but straightforward: 



'RSB 



n 

exp 



T^Va{p) 



ex^{-H[ipa{p)]} = 

d''pTT\n{p''6ab + rT, 



ab 



(3.9) 



The matrix under the logarithm in the above equation contains (k — l) eigen- 
values: 



X.=p' + r(l-^^^^l] (3.10) 
V uk- g J 

one eigenvalue: 

o / -u/c — 3q n\ 2uk r, 

X,=P^ + t[1- -^-^^l - r- 3.11 

y uk — g ) uk — g 

and in — k) eigenvalues: 

Substituting these eigenvalues into eq. (|3^ ), after simple algebra in the limit 
n one eventually obtains the following result: 
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Thus we see that in the region r ^ ^2/(4-d) ^^iq factor kg/{uk — g) in the 
exponential of the above equation is small compared to the leading term 
kT^'^~^^/'^ / (uk — g) given by the saddle-point solution, eq.( ^.l(j| ). 



4 Discussion 



It is interesting to note that non-analytic instanton contribution of the form 
given by eq. ( |2.21| ) can be easily "derived" basing on qualitative physical argu- 



ments. Let us again consider the random Hamiltonian (LI) at temperatures 
above Tc (r > 0), and let us estimate the contribution to the free energy 
coming from rare "ferromagnetic islands" where 6t{x) > t. In the mean- 
field regime at finite values of r the typical smallest (most probable) size of 
such island is Rc ~ r~^/^. Therefore, according to the probability distribu- 
tion, eq. (|1.2|) , in the limit of weak disorder {u —>■ 0) the contribution of the 
islands to the free energy with the exponential accuracy can be estimated by 
their probability: 

exp (-^^r(-^)/^) (4.1) 
which (up to undefined (const) factor) coincides with the result ( p. 211) . 



The above qualitative consideration seems rather valuable because it pro- 
vides good physical support for more exact but slightly formal and somewhat 
mysterious vector replica symmetry breaking scheme considered in Section 
2. 

Of course, exponentially small contributions to the free energy (as well as 
to others thermodynamical functions) of the type ( p.21|) are not so important 



for thermodynamical properties of the random ferromagnet in the considered 
paramagnetic temperature region. Nevertheless, the fact of their existence 
seems very interesting for two reasons. 

First of all, it tells that even in the mean-field regime the free energy 
of the random ferromagnet must be non-analytic function of the parameter 
which describes the strength of disorder m — >■ 0, which is interesting in itself. 
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Second, it indicates on the importance of non-linear excitations which in 
terms of the present rephca field theoretical approach are described by the lo- 
calized instanton-like solutions of the stationary equations. In the considered 
mean-field region away from these excitations provide only exponentially 
small corrections. However in the close vicinity of the critical point the pres- 
ence of instantons (which is ignored in the standard renormalization-group 
approach), and their interactions with the critical fluctuations may produce 
dramatic effect on the critical properties of the phase transition. It is wor- 
thing to note that although in the scaling regime (at T = Tc) the situation 
looks very different from that considered in this paper, the corresponding 
stationary equations ( p..8|) (with r = 0) also have non-linear instanton-like 
solutions with the RSB structure given by eq.( |2.5| ). One can easily check that 
in the dimension D = 4 these solutions can be found explicitly M: 



where the size of the instanton R appears to be the zero mode (the energy 
of the instanton does not depend on R). In dimensions below but close to 
four (at e = (4 — D) -C 1) the field configuration given by eq. (|4.2| ) can be 
considered as the approximate solution which contains the parameter R as 
the soft mode, since the energy of the instanton, eq. ( |2.71) , depends on R very 
weakly: 

E{k) = ^SnR-'-^ (4.3) 
3 uk — g 

(here So is the square of the unite D-dimensional sphere). 

At present it is not qute clear how all these non-linear instanton exci- 
tations could be incorporated into the self-consistent theory of the critical 
fluctuations. Keeping in mind that the degrees of freedom of this type ex- 
plicitly break the replica symmetry, a kind of "heuristic" renormalization 
group approach has been proposed 0, in which it was assumed that due to 
interactions of the fluctuations with this type of non-perturbative excitations 
the replica symmetry in the effective matrix, describing non-linear interac- 
tions of the fluctuating fields, is spontaneously broken. This resulted in the 
the instability of previously known fixed points and remarkable "runaway" 
behaviour of the renormalization group flows (which e.g. may indicate on 
the onset of a kind of the glass-like phase in a narrow temperature interval 
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around Tc). I hope the study described in the present paper would stimulate 
further much deeper investigation of the physics of the phase transition in 
random ferromagnets. 
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